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3.  ENERGY  ASPECTS  OF  RADAR  SCANNING  OF 
LIMITED  AREAS  OF  THE  EARTH'S  SURFACE 

Prof.  B.F.  Vysotskiy,  Doctor  of  Technical  Sciences, 
B.A.  Voynich,  Candidate  of  Technical  Sciences 


The  article  examinee  the  oharaoterietioe  of  energy  calcula¬ 
tions  for  airborne  radar  installations  operating  in  a  mode  of  lin¬ 
ear  scanning  of  the  earth's  surface .  It  considers  the  oases  of 
forward  scanning  and  lateral  scanning  in  relation  to  the  direction 
of  light.  The  width  of  the  scanning  strip  is  kept  constant  for 
different  observation  distances . 

Graphs  showing  the  signal-to-noise  ratio  as  a  function  of  the 
sighting  angle  of  the  scanned  areas  are  given. 


In  the  design  of  radar  sets  for  scanning  the  earth’s  surface 
(for  example,  sets  used  for  radar  cartography),  there  occur  cases 
which  require  scanning  of  the  surface  by  linear  plotting  of  its 
radar  image  with  a  fixed  constant  strip  width  bL.  In  general,  two 
types  of  scanning  may  be  distinguished: 

a)  scanning  in  which  the  scanning  ray  moves  in  the  direction 
of  flight  (forward  scanning); 

b)  scanning  in  which  the  scanning  beam  moves  sideways  from 
the  direction  of  flight. 

This  is  illustrated  in  Pig.  1,  which  shows  the  two  types  of 
linear  scanning  projected  onto  the  scanned  surface.  It  is  readily 
seen  that  in  linear  scanning  the  variation  of  echo-signal  power 
as  a  function  of  slant  distance  to  the  target  is  different  from 
the  variation  in  ordinary  cases. 

Using  the  fundamental  radar  equation 


we  find  that  the  echo-signal  power  for  a  fixed  direction  of  the 
antenna  in  a  horizontal  plane  is  determined  as  follows: 

liW'  (1) 
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where  ay  Is  the  width  of  the  directional  diagram  in  a  vertical 
plane. 

In  the  scanning  method  under  consideration  the  strip  width 
is  A L  ■  const.  Consequently,  as  the  distance  of  the  strip  from  the 
radar  set  increases,  its  complete  illumination  requires  a  narrower 
beam  than  the  illumination  of  a  nearby  strip.  However,  as  the  beam 
is  narrowed  (as  ay  is  reduced),  the  amplification  factor  of  the 

antenna  Increases.  The  Increase  in  amplification  factor  can  be 
compensated  by  a  decrease  in  the  signal  power  as  J?n  increases.  In 

other  words,  oy  is  a  function  of  Rn  which  can  Increase  as  J?n  in¬ 
creases. 


Pig.  1.  Patterns  for  linear  scanning  of  a  surface:  a)  scanning  in 
which  the  scanning  beam  moves  in  the  direction  of  flight  (forward 
scanning);  b)  scanning  in  which  the  scanning  beam  moves  sideways 
from  the  direction  of  flight.  The  arrow  indicates  the  direction 
of  flight,  and  the  dashed  line  indicates  the  trace  of  the  flight 
trajectory.  is  the  projection  onto  the  scanned  surface  of  the 

slant  distance  J?n  from  the  radar  set  to  the  edge  of  the  strip.  1) 

flgl;  2)  flg2;  3)  J?g3;  *0  radar  set. 


We  may  assume  that  both  the  flight  altitude  B  and  the  strip 
width  A L  are  specified  in  advance.  For  convenience  in  the  subse¬ 
quent  calculations,  we  shall  express  Rn  in  terms  of  H  and  the 

sighting  angle  0  of  the  strip: 

H  cosec  6  ( 2 ) 

and  we  also  express  ay  as  a„=a„(0)»  then 

p  ~ _ L _ _  (3) 

//«u*(G)c<ucc4fl  ’ 

In  calculations  involved  in  the  design  of  radar  sets,  it  is 
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more  convenient  to  use  not  the  absolute  value  of  the  echo  signal 
but  its  ratio  to  the  threshold  signal,  l.e., 


1 


q-Jjl _ 

x  VPm  V//*a;{6)co*ec<0 


1_ 


(4) 


The  purpose  of  the  present  study  is  to  explain  this  function.  In 
Fig.  2  and  in  the  text  we  have  adopted  the  following  notation:  H 
is  the  flight  altitude;  A£  is  the  width  of  the  scanning  strio  on 
the  earth's  surface;  J?g^,  J?g2,  ^3  are  the  horizontal  (earth)  dis¬ 
tances  to  the  first,  second,  and  third  scanning  strips,  respec¬ 
tively;  c*vl,  ay2,  ay3  are  the  widths  of  the  directional  diagram 

in  the  vertical  plane  which  are  required  for  illumination  of  the 
first,  second,  and  third  strips,  respectively;  yj,  Ya,  y*  are  the 
angles  between  the  antenna  beam  and  the  vertical  in  the  observa¬ 
tion  of  the  corresponding  areas  of  the  earth's  surface;  0  is  the 
sighting  angle  of  the  scanned  area  of  the  earth's  surface,  l.e., 
the  angle  between  the  antenna  beam  and  a  plane  parallel  to  the 
earth's  surface;  /?n  is  the  slant  distance  to  the  target;  Dy  is  the 

aperture  of  the  antenna  in  the  vertical  plane. 


PMC 


Fig.  2.  Diagram  showing  the  relative  positions 
of  the  radar  set  and  the  scanning  strip.  1) 
Radar  set;  2)  Fg2;  3)  /?g3;  4)  oyl;  5)  ay2J  6) 

°v3* 


Q  as  a  Function  of  the  $!gh+lng  Angle  0  of  the  Object  In  Forward 
Scann I ng 

We  shall  assume  that  the  radar  set  is  at  an  altitude  B  above 
the  surface  of  the  earth  (Fig.  2).  Then  in  order  to  observe  areas 
of  the  earth's  surface  which  have  a  width  of  AL  and  are  situated 
at  different  distances  from  the  radar  set,  we  can  use  antennss 
with  different  diagram  widths  in  the  vertical  plane. 

In  the  notation  we  have  adopted: 

Ya;  (5) 
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Rti  •)  ^ 
H 


~jj  tg  u»i ;  a*i  “  a  rctg  —  . 
■  tg  (Ya  -r  Ya J.-  0.3  “  a  rcig 


&a~-A— tKtVi+Ort)  Y» t  °»j arctg . 


(6) 

(7) 

(8) 


Prom  Formulas  (5)  and  (7)  it  follows  that 


=arctg 


f/tgva  -\-U 

ft 


-  Ya  “  a  rctg  ^tg  Vs  v  ~  j  ~  Ya-  (9) 


By  analogy  with  this,  from  Formula  (8)  we  obtain 


a,3  “  a  ret g 


H 


-  Yj—arctg  (tg  Yj  +  -  Yj.  ( 10 ) 


Let  us  now  write  the  expression  for  Q  for  forward  scanning 
in  the  case  of  observation  of  isotropically  scattering  point  ob¬ 
jects  resolvable  in  angle  and  range. 


We  assume 


V~Nrxli,  (U) 

(where  N.  is  the  number  of  object-echo  pulses  which  is  character- 

istic  for  the  ordinary  type  of  receivers  which  have  cathode-ray 
tube  radar  scopes  with  afterglow. 

When  the  object  sighted  is  situated  in  the  first  scanned 
area,  i.e.,  when  y  ■  Yi  *  0,  the  expression  for  Q 1  can  be  rewrit¬ 
ten  as 


since  — - .  If  we  write  hL/H  «  a,  the  formula  becomes 

COS  Hj!  ’ 


Qi 


A  cos-Uarclgw) 
ft*  (arctg<j)3 


(12) 


When  the  object  sighted  is  in  the  second  area,  i.e.,  when 
Y  j4  0,  the  formula  for  Q  is  obtained  by  substituting  into  (4)  the 
value  of  o  2  from  (9)  and  the  value  of  J?  for  the  case  y  /  0. 


h  _ h _ 

C0i  (Ya  +  a,a)  cos  [ys  +  arct"  (tg  va  *i-  a)  -  y2] 

^ _ H _ . 

cos  [arc!g  (tg  y-  v  a)l  ’ 

Q  ^  A  cos-»tarctg(tg  y2  a)] 

1  H*  [arctg(cg  yj  +  — Ysl2 


(13) 
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When  the  object  sighted  is  situated  in  the  third  area,  the 
expression  for  Q  can  be  obtained  from  (13)  by  a  simple  substitu¬ 
tion  of  the  angle  y*  for  the  angle  y*.  Prom  this  it  follows  that 
the  amount  by  which  the  assumed  power  exceeds  the  threshold  value, 
expressed  as  a  function  of  the  angle,  can  be  written  in  general 
form  as 


A_  con*  [arctg  (tg  y  «)[ 
tf*  (arctg  (tg  y  -:•«)  —  vlJ  ’ 


(14) 


If  in  (14)  we  replace  the  angle  y  by  the  angle  0  ■  90°  —  y, 
we  obtain  an  expression  for  Q  in  terms  of  the  sighting  angle  0  of 
the  object: 

Q  _  _/l_> _ cos<  [arctg  (Cts  C a)] _ 

1  H*  *  [a rctg  (ctg  0  a)  —  C0°  +  0J-  (15) 


The  variation  of  (15)  for  parameter  values  of  a  ■  10  and  a  ■ 
■0.25  are  shown  in  Pigs.  3  and  4. 

These  graphs  also  show  the  values  of  the  antenna  aperture  in 
the  vertical  plane,  Dyt  in  relative  units. 


•»  w  w  w  i 

Fig.  3.  Q  as  a  function  of  0,  the  sighting  angle  of  the  object, 
when  the  ratio  of  the  strip  width  A L  to  the  flight  altitude  H  is 

A L/H  m  a  -  10.  Curve  Qi  corresponds  to  forward  scanning,  V  N.  . 

-h  v 

Curve  Q 2  corresponds  to  lateral  scanning,  V  ^  .  Curve  Q s  cor¬ 

responds  to  lateral  scanning,  V  ^  N~. 1 .  The  target  is  a  point.  The 
diagram  is  symmetrical.  1)  Q  (in  relative  units);  2)  D y. 
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Pig.  4.  Q  as  a  function  of  6,  the  sighting  angle  of  the  object, 
when  the  ratio  of  the  strip  width  AL  to  the  flight  altitude  H  is 

_ it 

A L/H  *  a  «  0.25.  Curve  $i  corresponds  to  forward  scanning,  V  'v  N., 

-h  v 

Curve  Q 2  corresponds  to  lateral  scanning,  V  j?  .  Curve  Qs  cor¬ 
responds  to  lateral  scanning,  V  *  N~.1 .  The  target  is  a  point;  the 
diagram  is  symmetrical.  1)  Q  (in  relative  units);  2)  Z?v> 


Q  as  a  Function  of  the  Sighting  Angle  0  when  the  Scanning  Is  Lat¬ 
eral  In  Relation  to  the  Direction  of  Flight 

For  a  radar  set  with  lateral  scanning  we  can  also  derive  a 
function  similar  to  (15)  if  we  take  into  account  certain  charac¬ 
teristics  of  the  operation  of  such  sets. 

It  is  known  that  in  lateral  scanning  (or  when  the  antenna  is 
directed  sideways)  the  number  of  pulses  arriving  from  different 
objects  does  not  remain  constant  but  is  dependent  on  the  slant 
distance  of  these  objects. 

Nearby  objects  are  illuminated  by  the  set  for  a  shorter  time 
than  distant  objects. 


It  can  be  seen  from  Pig.  5  that  object  1  is  illuminated  by 
the  set  f o  a  length  of  time  t\ : 


*x 


A?ii  |fir 
vuyi 


and  object  2  is  illuminated  for 


•'nyi 


(16) 
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(17) 


where  ag  is  the  width  of  the  antenna  diagram  in  a  horizontal 
plane;  rput  is  the  ground  speed  of  the  aircraft. 

It  follows  from  Formulas  (16)  and  (17)  that  >  t\t  since 
*n2  >  *nl‘ 

Consequently,  the  number  of  pulses  from  the  object  depends 
on  the  slant  distance  to  the  object  and  is  expressed  by  the  for¬ 
mula 


Runr  p 

*  *  i% 

^nyt 


(18) 


where  F  .  is  the  repetition  frequency  of  the  radar  set. 

Therefore  the  value  of  the  visibility  coefficient  (11)  for  a 
radar  set  with  lateral  scanning  will  not  remain  constant  but  will 
depend  on  the  distance. 

The  formula  for  the  visibility  coefficient  may  be  written, 
on  the  basis  of  (18)  and  (11),  as 


a  i 

/>;  * 


(19) 


where  4i  Is  a  constant. 


Fig.  5.  Diagram  of  lateral  sighting  of  objects 
1  and  2  situated  on  the  surface  at  different 
distances  from  the  radar  set.  1)  flnl;  2)  J?n2; 

3)  radar  set;  4)  vput;  5)  otg;  6)  ay. 


In  order  to  obtain  an  expression  for  Q  for  a  radar  set  with 
lateral  scanning,  we  substitute  (19)  into  (4). 
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(20) 


Qi—A 


A  » » 


«•  —  «r^/a 

When  y  *  Yi  ■  0,  by  analogy  with  (12),  we  obtain 

n  _  .  _ 1 _  co«T;aaBi  __  A  ■/ toa7  ( a  fctg~o~) 

(arcii'j):*  //*  7/7/i’  (a  retort)2 

When  y  /  0,  the  expression  for  Q  Is  written  in  the  form 

O  -4  / cos2' f arcis  ( i.j  y  +  'a )~ 

3  [arctj(t«  y +  «)  — YlJ 


(21) 


(22) 


After  substituting  the  angle  0  ■  90°  -  y  for  the  angle  y  In  ( 22 ), 
we  obtain 


/  _ •/ cds7  [nrct;.  (i'iC'  H  '~<t 

[afcliT(7t-  C  i-  a)  «)’  6;  " 


*  »i 


(23) 


The  graphs  plotted  on  the  basis  of  Formula  (23)  for  parameter 
values  of  a  ■  10  and  a  *  0.25  are  shown  in  Figs.  3  and  4,  together 
with  the  graphs  of  (15).  The  calculations  using  Formulas  (15)  and 
(23),  necessary  for  the  plotting  of  these  graphs,  were  carried  out 
by  Engineer  D.V.  Terekhina. 


The  foregoing  calculations  show  that  in  linear  scanning  the 
scanning  of  distant  strips  is  most  advantageous  from  the  stand¬ 
point  of  energy.  The  strip  spacing  is  limited  by  the  design  diffi¬ 
culties  of  constructing  an  antenna  with  a  large  aperture  in  the 
vertical  plane. 

The  graphs  shown  in  Figs.  3  and  4  make  clear  the  energy  ad¬ 
vantages  of  lateral  scanning. 

Analogous  calculations  of  the  function  Q  «  f(Q)  for  cases  of 
targets  with  nonzero  length,  for  the  case  in  which  V  'v  N~. 1 ,  and 

also  for  cases  in  which  the  directional  diagram  used  is  such  that 

the  image  of  the  strip  is  uniformly  illuminated  have  shown  that 

the  laws  illustrated  in  Figs.  3  and  4  remain  valid.  As  an  example, 

Figs.  3  and  4  show  the  variation  of  Qj  for  the  case  in  which  V  * 

^  N~.1 . 

V 


Calculations  which  take  account  of  the  curvature  of  the  earth 
and  attenuation  in  the  atmosphere  were  not  carried  out.  However, 
it  may  be  assumed  that  the  effect  of  taking  account  of  the  curva¬ 
ture  of  the  earth  will  be  to  increase  Q  with  decreasing  0  (within 
the  limits  of  straight-line  visibility,  of  course),  while  taking 
account  of  attenuation  in  the  atmosphere  will  hinder  the  increase 
of  Q  with  decreasing  0. 


-  8  - 


Manu¬ 

script 

Page 

No. 

1 

1 

1 

2 

3 

6 


Transliterated  Symbols 

np  ■  pr  ■  priyemnik  ■  receiver 

h  ■  n  ■  naklonnaya  ■  slant 
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r  =  g  ■  gorizontal’naya  ■  horizontal 

in  ■  sh  ■  shum  ■  noise 

nyT  ■  put  ■  putevaya  *  ground 
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4.  SIMULATION  OF  RANDOM  NUMBERS  ON  AN  ELECTRONIC  DIGITAL  COMPUTER 

V.A.  Likharev,  Candidate  of  Technical  Sciences, 

Engineer  L.V.  Dobrolyubov,  Engineer  N.A.  Kobzev 


The  article  describes  a  procedure  for  obtaining  a  random  se¬ 
quence  of  numbers  with  a  given  distribution  law  by  using  uniformly 
distributed  random  numbers.  It  gives  programs  for  generating  the 
latter  on  the  BESM-2M  high-speed  electronic  computer.  As  examples , 
it  describes  the  methods  for  obtaining  uniform  normal ,  exponential , 
Rayleigh,  and  generalized  Rayleigh  laws ,  as  well  as  the  results 
obtained  by  checking  the  correlation  of  the  uniform  distribution , 
the  agreement  of  the  normal  distribution  with  a  given  distribution, 
and  the  estimation  of  numerical  characteristics  by  the  method  of 
confidence  intervals . 


With  the  advent  of  electronic  computers,  the  computation 
methods  grouped  under  the  general  name  of  Monte  Carlo  methods  are 
becoming  more  and  more  popular.  Two  tendencies  have  been  observed 
in  the  application  of  these  methods.  In  the  first  type  of  applica¬ 
tion,  the  abstract  mathematical  problem  of  solving  a  complicated 
equation  is  reduced  to  the  investigation  of  an  analogous  random 
process  described  by  the  same  relation.  For  example,  this  may  in¬ 
volve  a  partial  differential  equation  of  the  elliptic  or  parabolic 
type,  an  integral  equation,  or  the  calculation  of  a  multiple  in¬ 
tegral.  In  the  second  type  of  application,  a  given  random  process 
is  investigated  by  constructing  a  properly  idealized  mathematical 
model  of  it.  In  both  types  of  problems  the  necessary  characteris¬ 
tics  of  the  random  process  under  investigation  are  obtained  sta¬ 
tistically,  and  the  required  accuracy  is  achieved  by  increasing 
the  number  of  trials . 

This  last  fact  explains  why  such  calculation  methods  could 
come  into  use  only  after  the  appearance  of  electronic  computers, 
since  only  such  machines  are  capable  of  carrying  out  the  random 
process  under  investigation  a  sufficiently  large  number  of  times. 

With  a  sampled-data  digital  computer  it  is  possible  to  study 
the  value  of  the  process  at  successive  instants  of  time.  For  this 
purpose,  it  is  necessary  tc  obtain  by  means  of  an  electronic  com¬ 
puter  a  random  sequence  of  numbers  which  obey  some  distribution 
law  or  some  class  of  distribution  laws. 

The  methods  for  obtaining  random  sequences  of  numbers  which 
obey  given  distribution  laws  have  been  most  fully  discussed  in 
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[1],  The  normal  distribution  is  of  great  Importance  in  the  solu¬ 
tion  of  radar  detection  problems.  Let  us  consider  the  method  of 
obtaining  normal  numbers  by  means  of  electronic  computers. 


Suppose  that  it  is  required  to  obtain  a  sequence  of  random 
numbers  is.}  which  have  a  normal  distribution  with  a  mathematical 

t 

expectation  m  -  0  and  a  variance  o2  -  1, 


/.(*)  = 


VTn  o 


"RT 

e 


(1) 


According  to  the  central  limit  theorem  of  the  theory  of  prob¬ 
ability,  the  sums  of  a  large  number  of  random  terms  will  have  a 
distribution  which  tends  asymptotically  to  the  normal  distribution 
if  certain  very  general  conditions  are  satisfied.  Therefore  an  ap¬ 
proximate  simulation  of  normally  distributed  random  numbers  may  be 
obtained  by  summing  the  numbers  of  an  initial  sequence.  It  is  con¬ 
venient  to  use  quasiuniform  pseudorandom  numbers,  obtained  by  a 
programing  method,  as  the  initial  sequence. 

A  program  for  obtaining  quasiuniform  random  numbers  is  de¬ 
scribed  in  [2], 


However,  the  interval  of  aperiodicity  in  a  sequence  of  num¬ 
bers  obtained  by  means  of  thi3  program  does  not  exceed  50,000, 
which  may  be  insufficient  for  the  solution  of  a  number  of  prob¬ 
lems.  In  addition,  this  program  consists  of  four  commands,  and 
this  leads  to  considerable  waste  of  machine  time  when  the  opera¬ 
tions  must  be  repeated  a  large  number  of  times.  For  this  reason, 
the  program  we  used  as  the  initial-sequence  program  was  one  pro¬ 
posed  by  A. I.  Sragovich,  a  scientific  staff  member  at  the  Comput¬ 
ing  Center  of  the  Academy  of  Sciences  of  the  USSR  (see  Appendix). 

A  characteristic  of  this  program  is  the  fact  that,  first  of 
all,  it  consists  of  three  commands  and,  secondly,  that  Its  inter¬ 
val  of  aperiodicity  is  more  than  3* 10*.  The  numbers  obtained  by 
means  of  this  program  are  uniformly  distributed  over  the  interval 
(-1,  +1). 


The  basic  requirement  for  such  a  sequence  of  numbers  is  that 
there  must  be  a  low  statistical  correlation  between  its  numbers. 
The  lack  of  correlation  was  verified  as  follows.  We  calculated  the 
correlation  moments  of  the  sequence  by  means  of  the  formula 


ft- 1 


—y  X  l**  “  I  r*+/  ~  *1 


.  _  »-l 

*-i  \  *-i  / 


2  ' 


(2) 


where  n  is  the  number  of  terms  of  the  sequence;  l  is  the  number  of 

ft 

—  |  _ 

calculated  points;  £  x *  ls  the  mathematical  expectation  of 

A-l 

the  statistical  series;  x ^  are  the  values  of  the  terms  of  the  sta¬ 


ll  - 


tistical  series.  The  values  of  the  correlation  moments  calculated 
according  to  Formula  (2)  are  shown  in  Table  1. 


TABLE  1 


l 

01 

l 

01 

0 

1,000000000 

6 

-0,029306802 

I 

0,02-1834369 

7 

-0,025279736 

2 

—0,037379679 

8 

-0,0044158218 

3 

—0,027280937 

9 

-0,025902692 

4 

'  —0,037787475 

10 

-0,016323434 

5 

—0,025829599 

The  check  showed  that  the  correlation  properties  of  the  uni¬ 
form  numbers  do  not  change  over  the  block  of  numbers  from  0  to 
3* 10‘  . 


In  order  to  obtain  normal  numbers  from  the  initial  uniform 
sequence,  we  form  the  normalized  sum 


A 


h = k  S  5/ 1 

i-i 


(3) 


where  k  is  a  proportionality  factor  introduced  in  order  to  normal¬ 
ize  the  variance  a 2  =  1;  are  the  uniformly  distributed  numbers; 

n  is  the  number  of  terms . 


It  can  be  shown  that  h  has  a  nearly  normal  distribution  even 
for  relatively  small  values  of  n.  In  the  present  case  n  ■  5. 

As  was  shown  in  [4],  in  order  to  improve  the  asymptotic  nor¬ 
mality  of  ht  it  is  possible  to  use  special  transformations.  For 
example,  the  quantity 


V~h--±-(3h-h>)  (1») 

20n 

will  have  a  distribution  sufficiently  close  to  the  normal  distri¬ 
bution  for  much  smaller  values  of  n  than  in  the  case  of  simple 
summation. 

A  standard  program  for  obtaining  normal  numbers,  based  on  a 
uniform  distribution  and  making  use  of  Formulas  (3)  and  (4),  is 
given  In  the  Appendix. 

The  t  criterion  [3]  was  used  to  determine  whether  the  normal 
numbers  were  uncorrelated.  In  the  present  caee  this  criterion  is 
expressible  by  the  formula 


t 


~ _ Qi—  v 

V\  -of 


n-2, 


(5) 


where  p*  is  the  first  correlation  moment  of  the  normal  numbers. 
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The  value  of  pi,  calculated  for  a  block  of  10*  numbers,  was 
found  to  be 

q,  =  0,0042959261. 

If  the  value  of  t  obtained  by  means  of  Formula  (5)  exceeds 
the  tabulated  value  corresponding  to  the  assumed  significance 
level,  then  the  assumption  that  the  random  numbers  in  their  en¬ 
tirety  are  uncorrelated  is  incorrect.  In  the  case  of  a  5%  signifi¬ 
cance  level  for  the  sequence  under  investigation 

t  ■  0.136 

and  the  probability  p  of  finding  t  >  0.136  is  greater  than  5% • 
Consequently,  this  sequence  of  normal  numbers  may  be  considered 
uncorrelated. 

In  order  to  estimate  the  agreement  between  the  experimental 
and  theoretical  distributions,  we  investigated  the  resulting  num¬ 
ber  sequence  by  means  of  Pearson’s  x2  criterion. 

The  number  sequence,  consisting  of  10s  normal  numbers,  was 
divided  into  10  segments.  In  each  segment  we  took  the  first  104 
normal  numbers  and  plotted  a  histogram,  calculated  x2»  and  calcu¬ 
lated  m  and  o2  for  these.  The  results  of  the  calculations  are 
e  e 

shown  in  Table  2. 

TABLE  2 


<0 

>£ 

m. 

•l 

PaC.XOil.AC line  SKCIICpil- 
MCIIT,  II  TCOpCTII'l. 
paenp. 

2 

1 

0,0009335 

0,31031210 

24,378 

cayiaiiiioc  3 

2 

— 0,0024420 

0,32413262 

25,99 

cayviiiuoc  3 

3 

-0.0029806 

0,31925377 

27,3115 

c.iyin.iuGC  3 

4 

-0,0053246 

0,32229100 

36,911 

c.iyiaiiuoc  3 

5 

-0,0002083 

0,32151874 

22,14 

c.iyiaiiuoc  3 

6 

0,00001837 

0,32286972 

24,675 

c.iyiaiiuoc  3 

7 

-0,0038727 

0,31776674 

31,435 

cayiaiiimc  3 

8 

-0,0161672 

0,31574532 

36,138 

c.iyiaiiuoc  3 

9 

0,0024483 

0,31978057 

44,93 

cyuiccTuciiuoc  4 

10 

0,0034752 

0.3198C569 

42,725 

cyiuecTiiciiiioc  4 

1)  Number  of  segment;  2)  difference  between  ex¬ 
perimental  and  theoretical  distributions;  3) 
random;  b)  substantial. 


Figure  1  shows  a  histogram  of  the  first  10 4  normal  numbers, 
together  with  a  graph  of  the  theoretical  distribution.  Here  v ^  is 

the  number  of  normal  numbers  in  the  ith  segment  of  the  histogram. 
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In  estimating  the  agreement  between  the  theoretical  and  experimen¬ 
tal  distributions,  we  used  a  formula  from  [3].  If 


JLfclL 


>3. 


(6) 


the  difference  may  be  considered  substantial;  if,  on  the  other 
hand. 


"72v 


<3. 


(7) 


where  v  ■  r  -  1  ■  19,  r  is  the  number  of  segments  in  the  histo¬ 
gram,  the  difference  may  be  considered  random,  i.e.,  in  our  case. 


^KpUTH'l 


<37. 


Fig.  1.  Histogram  of  the  normal  distribution 
law:  v.  is  the  number  of  normal  numbers  in  the 

tth  segment  of  the  histogram.  1)  Theoretical 
distribution. 


As  can  be  seen  from  Table  2,  a  substantial  difference  on  the 
basis  of  the  x2  criterion  is  obtained  at  the  end  of  a  normally 
distributed  number  sequence  consisting  of  10®  numbers,  which  lim¬ 
its  the  range  of  numbers  used  to  9*10®  numbers. 

In  order  to  estimate  the  agreement  between  the  experimental 
and  theoretical  numerical  characteristics  of  the  normal  distribu¬ 
tion,  we  used  the  method  of  confidence  intervals.  Selecting  a  sig¬ 
nificance  level  of  If,  we  obtain  the  following  formulas  for  the 
boundaries  of  the  confidence  intervals  (see  [5]): 

tnt  -  -  2,58*„  <  ///,  <  ///r  f  2.583„,  (  8 ) 


where 


-  1H  - 


In  the  present  case,  for  «f 


-^_,/rtT~0  and  n  *  10 4 ,  we  have 


~0,0147</w,<  0,0147. 

The  confidence  interval  for  a  o^  with  the  same  level  of  reliabil- 

e 

lty  is  defined  as 

ot  -  2,58«.  + 2, 58j„  <9) 

where 


n 


In  the  present  case,  for  aT 


and  n  ■  10**, 


we  have  0.5583  £  o  <_ 


<  0.5817.  A  check  over  all  10  segments  showed  that  the  mathematical 
expectation  m  and  the  variance  0*  lie  within  the  99%  confidence  in¬ 
terval. 


The  normal  distribution  law  may  be  used  to  obtain  other  dis¬ 
tribution  laws  which  are  important  in  the  solution  of  radar  detec¬ 
tion  problems .  The  Appendix  gives  programs  for  obtaining  exponen¬ 
tial,  Rayleigh,  and  generalized  Rayleigh  distributions.  Each  of 
these  three  programs  includes  the  conversion  to  a  standard  program 
of  normal  numbers.  The  exponential  distribution  is  obtained  from 
the  normal  distribution  by  the  following  formula: 


TSZ  ^1  ho{>m 


The  formula  for  the  Rayleigh  distribution  is 

^I1**  E=  ^  NOPM 


(10) 

(11) 


The  generalized  Rayleigh  distribution  is  obtained  from  the 
normal  distribution  by  means  of  the  formula 

(12) 


where  a  is  a  constant. 

In  order  to  verify  the  quality  of  the  resulting  normal  se¬ 
quence  of  numbers,  we  solved  the  problem  of  detecting  a  nonfluc¬ 
tuating  target.  For  a  packet  of  10  pulses  and  a  false-alarm  prob¬ 
ability  of  F  *  10"6  we  obtained  three  values  of  the  probability 
of  correct  detection,  D>  corresponding  to  three  different  signal- 
to-nolse  ratios.  The  results  of  the  calculations  are  shown  in  Fig. 
2.  This  figure  also  shows  the  curve  obtained  by  calculation  in 
[6].  > 
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Pig.  ?.  Probability  of  detection  of  a  packet  of  nonfluctuating 

signals.  The  curve  was  calculated  for  the  parameters:  n  *  10,  F  * 

■  10“*;  A)  points  calculated  on  the  BESM-2M .  1)  U  ,  db. 

s 


The  standard  program  for  normal  numbers  makes  it  possible  to 
obtain  uncorrelated  normal  numbers  with  a  variance  a2  *  1  and  a 
mathematical  expectation  m  ■  0  and  may  be  recommended  for  use  with 
the  BESM-2M  high-speed  electronic  computer. 

The  programs  proposed  in  the  present  study  for  obtaining  ex¬ 
ponential,  Rayleigh,  and  generalized  Rayleigh  distributions  make 
it  possible  to  solve  a  number  of  radar  detection  problems. 

APPENDIX 

1.  Standard  program  for  obtaining  normal  numbers. 

Access  to  the  program  is  provided  by  V.M.  Kurochkin's  com¬ 
piler  program  by  the  command 

77  0000  0151  3777. 

The  result  -  the  normal  number  -  is  found  in  cell  0002.  This 
standard  program  can  operate  only  with  the  compiler  program. 

2.  Standard  program  for  obtaining  exponential  numbers. 

Access  to  the  program  is  provided  by  the  command 

34  0000  0027  0020. 

The  program  occupies  cells  0020  to  0030. 

The  quantity  y  -  q  +  1,  where  q 2  =  a|/cr*h  is  transferred  to 
cell  0003;  o*  i c  the  variance  of  the  signal  and  cr|h  is  the  vari¬ 
ance  of  the  noise. 
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The  result  is  found  in  cell  0004.  When  the  quantity  y  =  1  is 
sent  to  cell  0003,  the  result  obtained  is  a  sequence  of  exponen¬ 
tial  numbers  which  imitate  the  noise  voltages  at  the  output  of  the 
square-law  detector.  When  a  quantity  y  >  1  is  sent  to  cell  0003, 
the  result  obtained  is  a  sequence  of  numbers  which  imitate  the 
voltages  of  the  signal  together  with  the  noise  at  the  output  of 
the  square-law  detector  for  independent  fluctuations  of  the  sig¬ 
nals  reflected  from  the  target. 

3.  Standard  program  for  obtaining  Rayleigh  numbers. 

Access  to  the  program  is  provided  by  the  command 

34  0000  0031  0020. 

The  program  occupies  cells  0020  to  0032.  This  standard  program 
differs  from  the  previous  one  o"ly  in  the  fact  that  it  imitates 
the  voltages  of  the  noise  (when  y  =  1)  or  of  the  signal  together 
with  the  noise  (when  y  >  1)  at  the  output  of  the  linear  detector 
for  independent  fluctuations  of  the  signals  reflected  from  the 
target. 

i'he  quantity  y  =  q  +  1  is  sent  to  cell  0003,  and  the  result 
is  found  in  cell  0004. 

4.  Standard  program  for  obtaining  generalized  Rayleigh  num¬ 
bers  . 


Access  to  the  program  is  provided  by  the  command 

34  0000  0034  0020. 

The  program  occupies  cells  0020  to  0036. 

The  quantity  a  =  Ks//?agh  sen^  to  ceH  0003.  The  result 

is  found  in  cell  0004.  When  the  quantity  a  -  0  is  sent  to  cell 
0003,  the  result  obtained  is  a  sequence  of  Rayleigh  numbers,  and 
when  a  quantity  a  >  0  is  sent  to  cell  0003,  the  result  obtained 
is  a  sequence  of  numbers  distributed  according  to  the  generalized 
Rayleigh  law  and  imitating  the  voltages  of  the  signal  together 
with  the  noise  at  the  output  of  the  linear  detector  for  reflection 
from  a  nonfluctuating  target. 

Remark.  In  the  last  three  programs,  exit  from  the  program 
takes  place  upon  the  command  which  follows  the  program  access  com¬ 
mand. 
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•Commands  2002,  2003,  2004  constitute  A. I.  Sra- 
govich's  uniform-number  generator. 

1)  Standard  program  for  obtaining  normal  numbers:  2) 
supplement  to  kl ;  3)  uniform-number  generator*;  4)  con¬ 
stants  for  uni form- number  generator. 
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1)  Standard  program  for  obtaining  generalized 
Rayleigh  numbers. 
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8.  EFFICIENCY  OF  COHERENT  RADAR  IN  THE  RECEPTION 
OF  PARTIALLY  COHERENT  PACKETS 

Engineer  I.M.  Miroshnichenko 


The  present  study  determines  the  loss  at  the  threshold  sig- 
nal-to-noise  ratio  of  a  coherent  radar  in  the  reception  of  a 
packet  signal  with  a  frequency  different  from  the  fixed  tuning 
frequencies  of  the  Doppler  channels  or  with  a  frequency  which  var¬ 
ies  at  a  constant  rate. 

The  article  gives  formulas  and  graphs  by  means  of  which  it  is 
possible  to  select  the  radar  parameters  rationally  and  also  to  de¬ 
termine  the  requirements  for  the  stability  of  the  apparatus t 
starting  from  given  losses  and  probabilities  of  errors  (misses  and 
false  alarms)  for  coherent  and  coherent-noncoherent  treatment . 


A  great  many  published  studies  have  been  devoted  to  the  co¬ 
herent  detection  of  a  radar  signal,  but  in  general  these  studies 
either  did  not  take  phase  fluctuations  into  consideration  or  else 
treated  them  in  ways  which  required  the  use  of  a  complicated 
mathematical  apparatus  [5]. 

The  treatment  of  cases  of  regular  variation  of  a  Doppler  fre¬ 
quency  in  the  present  study  has  made  it  possible  to  compare  the 
efficiencies  of  coherent  and  noncoherent  reception  of  partially 
coherent  packets  in  a  fairly  simple  manner  and  to  formulate  the 
requirements  for  the  stability  of  the  phase-determining  elements 
of  the  apparatur. . 

The  use  of  coherent  accumulation  of  packet  pulses  in  ground- 
based  and  airborne  radar  sets  makes  it  possible,  as  is  known,  to 
detect  targets  at  much  greater  distances.  The  coherent  accumula¬ 
tion  is  carried  out  by  means  of  a  comb  filter  which  is  placed  in 
front  of  the  detector  and  whose  construction  (in  the  optimum  vari¬ 
ant)  involves  equal-weight  summation  of  the  packet  pulses  by  means 
of  a  set  of  delay  lines  [1],  In  practice  a  quasioptimal  accumula¬ 
tor  is  used  on  one  delay  line,  and  the  losses  are  small,  but  they 
increase  as  the  number  of  packet  pulses  increases. 

The  advantage  of  a  threshold  signal-to-noise  ratio  in  a  co¬ 
herent  optimal  accumulator  (which  sums  all  the  pulses  of  a  packet 
with  equal  weights)  increases  with  the  number  of  pulses.  From  the 
viewpoint  of  spectral  representations,  this  means  that  the  pass- 
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band  of  the  teeth  of  the  comb  filter  correspondingly  decreases. 

However,  as  the  number  of  pulses  in  a  packet  increases,  vari¬ 
ous  factors  which  lead  to  a  gradual  loss  of  coherence  come  into 
play. 


The  coherent  treatment  of  a  pulse  packet  is  optimal  for  a  co¬ 
herent  packet  against  a  background  of  Gaussian  white  noise,  i.e., 
a  packet  of  pulses  whose  initial  phase  varies  linearly  as  a  func¬ 
tion  of  the  number  of  pulses . 

The  derivative  of  the  phase  with  respect  to  time  is  the  Dop¬ 
pler  frequency  corresponding  to  the  radial  component  of  the  target 
velocity  with  respect  to  the  radar  set.  Since  the  target  velocity 
is  unknown,  a  coherent  system  consists  of  a  set  of  channels  M 
which  accumulate  pulses  for  fixed  Doppler  frequencies.  For  fre¬ 
quencies  different  from  these  fixed  values,  coherent  treatment 
will  result  in  some  loss  of  the  signal  energy.  If  the  velocity 
vectors  of  the  target  and  the  radar  set  do  not  lie  on  the  same 
straight  line,  the  frequency  of  the  reflected  signal  will  vary 
with  time,  and  the  relation  between  the  initial  phase  of  a  pulse 
and  its  number  will  be  nonlinear.  Consequently,  losses  in  a  co¬ 
herent  system  also  arise  as  a  result  of  failing  to  take  account  of 
relative  accelerations  of  the  target  motion. 

The  destabilizing  factors  disturbing  the  coherence  of  the  re¬ 
flected  packet  of  pulses  also  include  the  "blinking"  of  the  target, 
i.e.,  amplitude  and  phase  fluctuation  of  the  signal  as  a  result  of 
vibration  and  change  of  orientation  of  the  target,  which  will  usu¬ 
ally  be  of  fairly  complicated  form,  and,  finally,  instability  of 
the  radar  set  itself.  All  of  this  results  in  a  spread  of  the  spec¬ 
tral  components  of  the  signal  and  complicates  the  statistical 
properties  of  the  signal-noise  mixture.  A  coherent  treatment  of 
such  a  partially  coherent  packet  is  no  longer  optimal  and  leads  to 
an  increase  in  the  threshold  signal-to-noise  ratio  as  the  number  x 
of  pulses  in  the  packet  increases  (if  we  continue  to  narrow  the 
teeth  of  the  comb  filter). 

Since  the  signal  energy  losses  are  found  to  be  quite  consid¬ 
erable,  especially  for  high  target  velocities,  we  must  make  a  care¬ 
ful  examination  of  the  efficiency  of  coherent  systems  and  work  out 
recommendations  for  the  selection  of  their  fundamental  parameters. 
The  results  obtained  will,  of  course,  also  apply  to  coherent  fil¬ 
ter  systems  [2], 

From  this  point  on  we  shall  assume  that  there  is  an  optimal 
system  for  detecting  a  coherent  packet  of  a  signal  with  an  unknown 
(random  and  equiprobable)  initial  phase  against  a  background  of 
Gaussian  white  noise  with  two  squaring  channels  [3],  [53.  ) 

We  assume  that  the  radar  set  emits  a  sequence  of  coherent 
pulses  (or  has  a  coherent  heterodyne  oscillator  phased  by  the 
transmitter)  and  scans  the  space  in  some  sector  with  a  narrow 
beam.  Assuming  that  the  target  moves  in  a  straight  line  with  a 
constant  velocity  W  with  respect  to  the  radar  set,  we  shall  denote 
the  distance  between  them  by  R(t)  and  the  angle  between  the  radius 
vector  from  the  target  to  the  radar  set  and  the  velocity  vector  W 
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by  6. 


As  the  target  moves  during  the  tiyie  it  passes  through  the  an¬ 
tenna  beam,  the  coherence  of  the  reflected  packet  of  pulses  is 
disturbed.  In  a  manner  analogous  to  what  was  done  in  [4],  we  shall 
expand  the  function  R(t)  during  the  period  of  observation  into  a 
power  series  in  t  -  to  (where  to  is  the  start  of  the  observation) 

H  (0 = R  (/0)  -i-  ('  -  O  R'  (Z0)  (/  O'  !<"  I O 


where 


R’  (/0)  —  W  cos  0; 

A’"(/o)  =  “Sii:=0. 

r\ 

Accordingly,  the  initial  phase  of  the  reflected  signal  will  have 
the  form 


?a(0-1£^cos  0>(/ sh.J 0-{/  -  /„)*. . . 

A  '  r\ 


(1) 


(where  terms  of  higher  than  second  order  in  the  above  series  are 
neglected).  The  second  term  is  obviously  the  parasitic  displace¬ 
ment  of  the  initial  phase  which  leads  to  loss  of  coherence.  Taking 
account  of  the  packet  structure  of  the  reflected  signal,  we  write 
the  expression  for  the  parasitic  phase  displacement  of  the  feth 
pulse  in  the  form 


where 


o~ 


2nt V'-T 

~TfT 


2 

-  sin2  G. 


The  coefficient  a  determines  the  displacement  of  the  initial 
phase  between  the  first  ( k  ■  0)  and  the  second  ( k  *  1)  pulses  of 
the  packet  and  depends  not  only  on  the  course  of  the  radar  set  and 
the  target,  the  velocities,  the  difference  between  their  altitudes, 
and  the  target  distance,  but  also  on  the  parameters  of  the  radar 
set,  X  and  The  parasitic  displacement  of  the  Doppler  frequency 

is  of  the  form 


The  value  of  the  coefficient  a  and  the  corresponding  rate  of 
change  of  the  Doppler  frequency,  dF^/dt,  for  a  few  particular 

cases  are  shown  in  Table  1  [8]. 


Let  us  consider  the  radial  motion  of  the  target  when  the  Dop¬ 
pler  frequency  is  different  from  the  fixed  frequency  values  to 
which  the  Doppler  channels  are  tuned.  In  this  case  we  may  assume 
that  the  kth  pulse  of  the  reflected  packet  has  an  initial-phase 
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displacement  cp^,  which  is  due  to  the  difference  between  the  fre¬ 
quencies  of  the  reflected  signal  and  the  nearest  channel.  For  a 
frequency  midway  between  the  frequencies  of  two  adjacent  channels, 
the  phase  value  cp^  will  be  maximum.  Let  us  consider  this  case. 

Let  the  range  of  Doppler  frequencies  De  A and  let  the  number  of 

channels  be  M.  Then  the  frequency  which  falls  midway  between  chan¬ 
nels  m  and  m  +  1  is 


Fu 


(2) 


where  m  is  the  number  of  the  channel,  m  “  0,  1,  2  ...  Mt  and  its 
difference  from  the  channel  frequencies  is 


K 


To  this  Doppler  frequency  will  correspond  a  displacement  of  the 
initial  phase  of  the  kth  pulse  which  is  not  compensated  in  the 
channel, 

(3) 

where 


x.w 

A W  is  the  range  of  target  velocities. 

Table  2  gives  the  values  of  b  for  a  few  particular  cases. 
(The  physical  meaning  of  the  coefficient  b  is  the  displacement  of 
the  initial  phase  from  pulse  to  pulse.) 

In  the  general  case  of  relative  motion  of  the  radar  set  and 
the  target,  the  phase  of  the  fcth  pulse  is  determined  on  the  basis 
of  (1)  as 


yu-bk  +  tk-.  (2|) 

In  order  to  simplify  the  picture,  let  us  subdivide  the  general 
case  into  two  simple  cases,  each  corresponding  to  one  term  of  this 
expression. 

The  first  case  [described  by  the  first  term  of  the  expression 
(4)]  is  a  radial  displacement  of  the  target,  where  the  Doppler 
frequency  lies  midway  between  two  channels  and  is  described  by  the 
expression  (2). 

The  second  case  (described  by  the  second  term)  is  a  nonradial 
motion  with  an  initial  frequency  value  matching  one  of  the  chan¬ 
nels,  e.g.,  channel  m. 

After  determining  the  losses  at  the  threshold  signal-to-noise 
ratio  for  each  of  these  cases  individually,  we  can  easily  turn  to 
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a  consideration  of  the  more  general  picture  of  relative  motion  of 
the  target  with  respect  to  the  radar  set. 

To  determine  the  loss  at  the  threshold  signal-to-noise  ratio 
in  coherent  treatment  of  partially  coherent  packets ,  we  can  use 
the  following  method.  A  packet  of  partially  coherent  pulses  with 
an  arbitrary  envelope  can  be  reduced  under  certain  conditions  to 
an  equivalent  coherent  packet  with  the  same  number  of  pulses  and 
a  rectangular  envelope,  with  some  effective  signal  level.  For  this 
reason  it  seems  possible  to  use  the  theory  of  optimal  treatment 
methods  worked  out  for  such  a  packet.  In  particular,  it  is  pos¬ 
sible  to  use  formulas  and  graphs  for  the  probabilities  of  correct 
detection  and  of  a  false  alarm.  For  this  purpose,  we  must  obvi¬ 
ously  have  an  equivalent  packet  at  the  output  of  the  apparatus  un¬ 
der  consideration  with  exactly  the  same  statistical  properties 
(multidimensional  distribution  law  or  characteristic  function)  as 
the  original  packet.  For  the  sake  of  simplicity,  we  shall  consider 
rectangular  pulse  packets  which  fluctuate  in  unison.  Let  us  find 
first  of  all  the  voltage  at  the  output  of  the  intermediate-fre¬ 
quency  amplifier,  regarding  it  as  an  optimal  filter  for  each  in¬ 
dividual  pulse.  Let 


the  signal  of  the  feth  pulse  of  a  pure  coherent  packet  with  a  ran¬ 
dom  and  equiprobable  phase  a  at  the  input  of  the  receiver. 

Let  us  write  the  expression  for  the  signal  of  the  fcth  pulse 
of  a  partially  coherent  packet  at  the  input  of  the  receiver  in 
the  form 


OV-T^-T"). 

where  we  use  to  denote  the  initial  phase,  which  we  shall  assume 


the  second  case. 


We  assume  that  the  condition 


“0 


where  t  is  the  duration  of  the  pulse;  u>o  is  the  (intermediate) 
filling  frequency;  n  is  an  integer,  is  satisfied. 

Such  a  mathematical  mocjel  presupposes  that  for  sufficiently 
short  pulses  the  Doppler  effect  will  be  manifested  only  in  the 
change  of  the  initial  phase  of  the  pulses. 


After  filtering,  these  signals  at  time  t  =  t  will  be  equal, 
respectively,  to 


u:ii~  — 


0 


„  cl: 

tut  - - -  - —  cos  a---  •  cos  a; 

-<'0  Gq 
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where  Go  is  the  density  of  the  energy  spectrum  of  the  noise;  E  is 
the  energy  of  an  individual  signal  pulse;  uo>  are  the  intermedi 

ate  and  high  circular  frequencies;  a  is  a  dimensionality  coeffi¬ 
cient  . 

After  the  intermediate-frequency  amplifier  the  signal  and  the 
noise  go  to  two  squarers  which  have  synchronized  detectors,  accu¬ 
mulators,  and  squarers,  after  which  the  voltages  of  the  two  chan¬ 
nels  are  added.  Nonlinear  operations  -  extraction  of  the  square 
root  and  In  lots)  -  are  performed  on  the  sum,  and  this  is  followed 
by  the  threshold  apparatus.  We  find  the  signals  at  the  output  of 
the  linear  part  of  the  system. 

After  the  synchronized  detectors  and  accumulators  in  the 
squaring  channels  the  signals  u2fc  and  are  eQual»  respectively 

to : 

\r  rf- 

Uy-~-  N  •  -co;u 
fie 

Go 


where  N  is  the  number  of  pulses  in  the  oacket;  u,  V  are  the  volt 
ages  in  the  first  and  second  squaring  channels; 


A’ 


A 

cos 

A  t 


«.*  /: 


After  squaring  the  voltage  in  each  of  the  squaring  channels, 
adding,  and  taking  the  square  root,  we  find  the  envelope  under 
which  the  operation  In  Io(x)  is  carried  out.  If  we  equate  the  en¬ 
velopes  of  the  oscillations  at  the  outputs  of  the  accumulators  of 
the  coherent  and  partially  coherent  packets,  the  distribution  laws 
for  these  two  packets  will  be  equal,  and  consequently,  so  will  the 
probabilities  of  error  after  optimal-system  treatment.  Therefore 
the  effective  signal  level  must  be  found  from  the  conditions  of 
equality  of  the  envelope. 

The  expressions  for  the  envelope  of  the  coherent  and  partial¬ 
ly  coherent  packets  against  a  noise  background  at  time  t  »  t  are 
written,  respectively,  in  the  forms  [6] 


Go 


//(t)  -j- 


rli  ..  i- 

•  A  s: . . j 
fi«  I 
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where  A(t )  t  B( t)  are  independent  Gaussian  functions  with  a  correla¬ 
tion  coefficient  pfi)  and  a  variance  a  (noise  voltage).  If  we  in¬ 
troduce  the  notation 


X~y..  *«««.».;  | 

•  0  0  (t01  i 

then  the  above  formulas  can  be  rewritten  as 


X,  —  \f  l«  [  </  .V  cu>  ct ] -  — \b-\-  V  qN  .sin  ct)-; 


.u-}/ -;-r </V cos (u \  </ Vsin(«-l-oi)|3, 


equating  xi  and  xi  and  passing  to  the  effective  signal-to-noise 
ratio  of  a  coherent  packet  qe%  we  obtain  the  following  conditions: 

.v 

|  </ArCOS  u  —  |  q .  V  cns(u-|  V); 

*  1 

A' 
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Eliminating  the  parameter  a  and  assuming  that  the  signal-to-noise 
ratios  are  threshold  values,  we  obtain  the  final  expression 


wo? 


V  I.Op* 


(5) 


where 


For  the  first  case  of  target  motion  under  consideration,  we  have 
h  *  ^kl* 

For  sufficiently  large  values  of  N  we  can  replace  the  summa¬ 
tion  sign  with  an  integral  sign.  After  some  fairly  simple  trans¬ 
formations,  we  obtain  the  expression 


bX 
sir.  -- 


0 


(6) 


If  the  initial  phase  of  the  signal  is  known  (a 

tain 


sin  l>X 


0) ,  we  ob- 
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manner,  we  can 


For  the  second  case  ^  «  <pfe.  In  an  analogous 
obtain 

Y1-  - -r  I  /  -f  rCWvW 

/V  J/  2.i 

where  C(aNz)  and  S(aN 2)  are  Fresnel  Integrals. 

For  a  signal  with  a  known  initial  phase  (a  ■ 
mula  and  the  succeeding  formulas,  the  second  term 
root  sign  must  be  set  equal  to  zero. 

The  curves  plotted  in  accordance  with  Formula  (6)  for  three 
values  of  b:  b\  ■  0.21;  b 2  ■  0.08  and  2? 3  ■  0.04  radians,  in  accord 
ance  with  Table  2,  are  shown  in  Fig.  1. 

1 

t* 
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Fig.  1.  Loss  coefficient  as  a 
function  of  the  number  of  pulses 
in  the  packet  when  the  target 
moves  radially. 


The  curves  plotted  in  accordance  with  Formula  (7)  for  three 
values  of  a:  a\  ■  l.TS’lO"*,  a%  ■  3*10"4  and  as  ■  ^*10”5  radians, 
in  accordance  with  Table  1,  are  shown  in  Fig.  2. 


(7) 

0),  in  this  for- 
under  the  square 


Fig.  2.  Loss  coefficient  as  a 
function  of  the  number  of  pulses 
in  the  packet  when  the  target 
moves  nonradially. 
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Now,  in  order  to  plot  the  desired  characteristics  of  Figs.  3 
and  we  can  use  the  graphs  given  in  [7],  which  give  the  relation 
between  the  probabilities  of  error  Ppr  and  P ^  ^  for  different 

threshold  signal-to-noise  ratios  in  an  optimal  detection  system 
for  a  packet  of  coherent  pulses  with  phases  fluctuating  in  unison, 
as  well  as  Kaplan's  formula,  as  corrected  in  [3],  for  an  optimal 
detection  system  for  packets  of  incoherent  pulses. 


Fig.  3.  Threshold  signal-to- 
noise  ratio  as  a  function  of 
the  number  of  pulses  in  the 
packet  when  the  target  moves 
radially.  1)  «por. 


Fig.  4.  Threshold  signal-to- 
noise  ratio  as  a  function  of 
the  number  of  pulses  in  the 
packet  when  the  targest  moves 
nonradially.  1)  <7  . 


It  is  shown  in  [3]  that  for  low  error  probabilities  there  ex¬ 
ists  an  approximate  relation  (with  an  error  of  0.5  db)  between  the 

P„  Pi  4. ,  and  <7 _ of  a  coherent  packet: 

p.r’l.t  ^por 


/ 


!u 


This  characteristic  is  represented  in  Figs.  3  and  *1  by  curve 
2.  For  comparison,  these  figures  also  show  the  characteristic 
(curve  1)  for  incoherent  treatment  of  the  packets  (plotted  accord- 
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ing  to  Kaplan’s  formula).  Curves  3,  **,  and  5  (Pigs.  3  and  4)  for 
coherent  treatment  of  partially  coherent  packets  are  plotted  on 
the  basis  of  Pigs.  1  and  2  and  Formula  (5),  with  the  curves  of 
Pig.  3  corresponding  to  the  curves  of  Pig.  1  and  the  curves  of 
fig.  S  corresponding  to  the  curves  of  Pig.  2.  In  general,  the 
curves  of  Pigs.  3  and  4  are  similar.  As  the  number  of  pulses  in 
the  packet  increases,  the  efficiency  of  coherent  treatment  rapidly 
decreases  if  we  do  not  take  account  of  the  relative  acceleration 
of  the  target  or  if  we  use  an  insufficient  number  of  Doppler  chan¬ 
nels. 


Physically  this  is  explained  by  the  fact  that  as  N  increases, 
the  parasitic  displacement  of  the  initial  phase  also  increases 
from  pulse  to  pulse.  The  accumulation  of  such  out-of-phase  pulses 
leads  to  inevitable  losses  in  signal  energy,  while  the  noise  power 
continues  to  increase  with  increasing  N.  For  a  linear  displacement 
of  the  initial  phase  (Pigs.  1,  3)  it  is  possible  to  compensate  for 
the  accumulated  pulses  almost  completely,  and  therefore  the  maxima 
of  the  curves  of  Pigs.  1  and  3  may  be  large. 

In  considering  the  second  case  of  target  displacement  we  as¬ 
sumed  that  the  initial  value  of  the  Doppler  frequency  (of  the  re-? 
ceived  packet)  corresponds  to  channel  m  and  that  its  variation 
causes  a  loss  of  signal  energy  in  this  same  channel.  However,  as 
the  duration  of  the  packet  increases,  the  maximum  value  of  the  ao-? 
cumulated  voltage  passes  from  channel  m  to  channel  m  +  1,  then  to 
channel  m  +  2,  and  so  on. 

The  loss  of  signal  energy  will  increase,  but  not  as  fast  as 
in  channel  m.  We  must  therefore  determine  what  is  the  loss  in 
channel  m  when  the  initial  value  of  the  Doppler  frequency  falls 
into  I  channel  [sic]. 

In  the  notation  of  Formula  (4)  the  initial  phase  of  the  feth 
pulse  in  channel  m  is 


<Pimi  =••  2  b.;m—ak'\ 

For  large  values  of  N 


After  some  fairly  simple  transformations,  we  can  reduce  this 
expression  to  the  form 


where  C  and  S  are  Fresnel  integrals . 
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This  makes  it  possible  to  determine  the  coefficient  y (  m) * 

and  consequently  also  the  threshold  signal-to-noise  ratio  for  each 
channel  of  a  coherent  system. 

For  example,  curve  6  (Fig.  4)  was  plotted  for  the  second 
channel  (m  =  1),  with  b  ■  0.04,  a  -  1.75*10“3  and  M  a  1000.  How¬ 
ever,  when  the  number  of  channels  is  of  the  order  of  hundreds  (and 
when  their  construction  is  optimum),  it  is  possible  to  use  curves 
3,  4,  and  5  of  Fig.  4,  neglecting  the  adjacent  channels. 

The  expression  (8)  makes  it  possible  to  construct  curve  7 
(Fig.  4),  which  gives  a  lower  bound  for  the  possible  losses,  for 
example,  at  a  *  1.75*10“3.  For  this  purpose  we  take  the  maximum 
value  of  the  Fresnel  integrals  (0.78),  as  well  as  the  maximum  val¬ 
ues  of  the  sine  and  cosine  (1).  Then  (l/y2)mln  =  ^2/8800,  i.e., 

the  losses  due  to  the  relative  acceleration  of  the  target  cannot 
be  less  than  a  fixed  value  even  if  there  are  infinitely  many  chan¬ 
nels.  These  losses  will  increase  as  N2 ,  and  consequently  for  any 
parameters  of  the  radar  set,  starting  from  some  value  of  N  depend¬ 
ing  on  the  value  of  the  relative  acceleration  of  the  target,  a  co¬ 
herent  system  will  become  less  efficient  than  an  incoherent  system. 

Using  the  graphs  plotted  in  Figs.  3  and  4,  we  can  make  a  ra¬ 
tional  choice  of  the  number  of  channels  and  the  number  of  delay 
lines  of  the  accumulator  (the  passbands  of  the  teeth  of  the  comb 
filter).  Let  us,  for  example,  choose  the  number  of  channels  in  ac¬ 
cordance  with  curve  5  (Pig.  3).  We  shall  assume  the  loss  for  a 
Doppler  frequency  lying  midway  between  the  channels  to  be  equal  to 
the  segment  6.  Since  this  loss  value  corresponds  to  a  number  N o  of 
pulses,  the  number  of  delay  lines  of  the  accelerator  must  also  be 
equal  to  N 0,  and  the  passband  of  the  teeth  of  the  comb  filter  [1] 
is 


2Si- 


_fjL 

2iV0  ‘ 


The  graph  of  <7por  for  the  fixed  Doppler  frequencies  corres¬ 
ponding  to  the  tuning  of  the  channels,  if  N  >_  N o,  will  be  deter¬ 
mined  by  the  straight  line  7.  For  other  values  of  the  Doppler  fre¬ 
quencies,  <7por  will  vary  within  the  limits  of  the  segment  6. 

The  threshold-signal  advantage  over  an  incoherent  accumulator 
will  decrease  as  N  increases  and  will  become  a  disadvantage. 

Because  of  these  losses  caused  by  the  gradual  violation  of 
the  coherence  of  the  reflected  signal,  it  is  of  considerable  in¬ 
terest  to  consider  a  mixed  type  of  accumulation,  in  which  the 
packet  is  divided  into  groups  of  coherently  accumulated  pulses, 
in  accordance  with  a  selected  coherence  interval,  and  the  pulses 
of  all  the  groups  are  then  accumulated  incoherently. 

If  we  take  the  coherence  interval  initially  to  be  No  ,  then  in 
Fig.  3  we  can  plot  graph  8  for  mixed  accumulation  of  the  signal, 
assuming  that  the  values  on  the  abscissa  axis  represent  not  fhe 
number  of  pulses  in  a  packet  but  the  number  of  pulse  groups  accu- 
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mulated  incoherently.  Graph  8  is  plotted  in  the  same  way  as  graph 
1,  according  to  the  corrected  Kaplan  formula  but  taking  account  of 

the  new  value  of  q  „  for  N  ■  1.  A  similar  construction  for  mixed 

Por 

treatment  can  also  be  carried  out  in  Pig.  4. 

It  is  clear  from  the  figure  that  whenl the  number  of  coherent 
pulse  groups  in  the  packet  is  small  (less  chan  5)»  the  mixed 
treatment  of  a  partially  coherent  packet  is  almost  as  efficient  as 
the  coherent  treatment  of  a  strictly  coherent  packet  with  the  same 
number  of  pulses.  In  the  general  case  -  for  any  N  -  the  threshold 
signal-to-noise  ratio  for  mixed  treatment  is  always  less  than  for 
incoherent  treatment  by  the  value  of  q  _  of  the  accumulated  co- 

herent  group  (determined  by  the  interval  of  coherence). 

The  foregoing  analysis  makes  it  possible  to  formulate  the  re¬ 
quirements  for  the  stability  of  the  phase-determining  elements  of 
a  radar  set:  the  transmitter,  the  local  and  coherent  heterodyne, 
and  the  delay  lines.  It  seems  possible  to  take  into  account  both 
the  absolute  value  of  the  frequency  difference  A /  (the  frequency 
drift)  and  the  drift  rate  A//At. 

The  value  of  A /  is  easy  to  find  from  Fig.  3,  given  the  allow¬ 
able  loss  in  the  threshold  ratio  q  .  For  this  it  is  necessary  to 

^  por 

determine  from  the  curves  the  required  kdop  for  a  given  N.  Then 
the  desired  requirement  is  written  in  the  form 

I  A/  U„  . 

•  U 

The  value  of  A //At  can  be  found  from  Fig.  4,  again  after 
specifying  the  loss  in  <7por  and  determining  for  the  required  N  the 

necessary  value  of  the  coefficient  adop*  The  desired  requirement 

is  of  the  form  |A//At|  <  adop/7T2’p. 

Thus,  in  this  article  we  have  considered  a  method  for  deter¬ 
mining  the  signal  threshold  power  loss  arising  in  a  coherent  sys¬ 
tem  both  because  of  the  limited  number  of  channels  and  because  of 
the  relative  acceleration  of  the  target  in  nonradial  motion.  The 
analysis  shows  that  even  though  the  advantage  of  a  coherent  system 
increases  as  the  number  of  pulses  in  a  packet  increases,  it  be¬ 
comes  more  and  more  difficult  to  attain  this  advantage  in  practice. 
It  is  necessary,  to  use  a  much  larger  number  of  channels.  Further¬ 
more,  at  the  present-day  velocities  of  targets  and  radar- carrying 
vehicles,  the  accumulation  necessarily  requires  the  introduction  of 
a  phase  correction  for  the  variation  of  Doppler  frequency  due  to 
the  relative  acceleration  of  the  target,  since  even  an  unlimited 
increase  in  the  number  of  channels  tuned  to  constant  Doppler  fre¬ 
quencies  will  not  make  possible  any  substantial  reduction  of  the 
signal  threshold  power. 

The  article  gives  formulas  and  graphs  which  can  be  used  for 
selecting  the  number  of  channels,  the  passbands  of  the  comb  filter 
and  the  other  parameters  of  the  radar  set  for  different  relative 
accelerations  of  the  target }  determining  the  requirements  for  the 
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stability  of  the  apparatus,  comparing  different  forms  of  treat¬ 
ment  -  incoherent,  coherent,  and  mixed  treatment  -  from  the  view¬ 
point  of  threshold  power,  and  selecting  the  interval  of  coherence 
for  given  probabilities  of  correct  detection  and  false  alarms. 
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Transliterated  Symbols 

a  =  d  =  dopplerovskaya  =  Doppler 
n  =  p  =  povtoreniye  =  repetition 
b  =  v  =  vertikal ’naya  =  vertical 
a  =  e  =  effektivnoye  =  effective 
nop  =  por  =  porogovoye  =  threshold 
np  ■  pr  =  pravil’noye  =  correct 
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ji.t  ■  l.t  ■  lozhnaya  trevoga  ■  false  alarm 


31 

31  n.p  ■  p.r  «  pravil'noye  ■  correct 

3I1  flon  ■  dop  *  dopustimyy  »  allowable 
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